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We show that a single I — 1 spin-parity J PC — 1 ++ a\ resonance can manifest itself as two 
separated mass peaks, one decaying into an S-wave pn system and the second decaying into a P- 
wave /o( 980 ) 7 t system, with a rapid increase of the phase difference between their amplitudes arising 
mainly from the structure of the diffractive production process. This study clarifies questions related 
to the mass, width, and decay rates of the a\ resonance raised by the recent high statistics data of 
the COMPASS collaboration on a\ production in ttN —»• ttttttN at high energies. 


PACS numbers: 12.40.Yx, 13.25.Jx 

New insight into the properties of light mesons is 
emerging from the unprecedented statistical precision of 
the COMPASS experiment at CERN where beams of 190 
GeV pions interact with nucleon targets [l|. These data 
are bound to enrich (if not challenge) our understand¬ 
ing of low-energy meson spectroscopy while, in addition, 
uncovering possible evidence for long-sought states of 
the strong interaction QCD potential beyond the quark- 
antiquark states of the standard model. 

We focus on the isospin 1 axial-vector resonance 
a% (reported in the Particle Data Compilation as 
a(l)(1260) i). Evidence is presented in the COMPASS 
data for a new narrow J PC = 1 ++ axial-vector state, 
strongly coupled to the 7t/o(980) system. This observa¬ 
tion of a peak in the two body 7r/o P-wave intensity at a 
mass of 1.42 GeV, combined with a phase motion close to 
180° with respect to other waves, appears at face-value 
to mean that a second axial-vector resonance is present, 
close in mass to the known broad ai(1260) that couples 
mainly to the np system [1]. While these three features, 
i.e. two peaks at different masses and a rapid phase varia¬ 
tion, are clearly present, there are reasons to be surprised, 
among which we mention: (i) The ai(1260) is a central 
member of the axial-vector nonet, which, together with 
the J p = 0 - , l - , and 0 + form the ground-state of the 
light quark-antiquark spectrum. A newcomer in the fam¬ 
ily would be difficult to accommodate, (ii) It is peculiar 
to have two J PC = 1 ++ three-pion states, with identical 
quantum numbers, close in mass (within a full width of 
each other), with orthogonal decay modes, without the 
presence of some new quantum number. The — K\ 
system led to decisive discoveries in fundamental physics; 
neutrino mixing is a spectacular current example. How¬ 
ever, in the a\ case, we see no candidate for a distinguish¬ 
ing quantum number. 

Our basic approach to high energy forward production 
of three pion states in pion-nucleon interactions is the 
Drell-Hiida-Deck mechanism [3]. This model has been 


studied extensively in the production of the J p = 1 + 
pn system Q, and here we extend the analysis to the 
J p = 1 + /o(980)7r system. An important difference is 
that whereas the pir system is in an orbital S-wave state, 
the /o(980)7r is in an orbital P-wave state. Since the 
two-body pir and fair systems are strongly interacting, 
we must modify the Deck mechanism with the proper 
final-state interactions due to the re-scattering of these 
systems. This is an inescapable physical consistency con¬ 
dition of the entire analysis. The unitary coupled chan¬ 
nel approach that we developed in Refs. m3 , should be 
an ideal way to show whether one resonance suffices or 
whether the COMPASS data do require two nearby res¬ 
onances with the same axial-vector quantum numbers in 
the three-pion system. 

In this Letter, we demonstrate that a single resonance 
suffices to explain the data and that the fair decay mode 
of the usual a\ is being observed for what appears to be 
the first time [2|. Our method can be used to determine 
new values for the mass and width of the a\, information 
important for lattice QCD and other calculations of the 
hadron spectrum. 
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FIG. 1: Deck production processes for (a) pn and (b) fon. 

Two-channel Deck amplitudes. We follow closely 
Refs. Q and We consider the Deck amplitudes 

for production of the quasi-two body systems np and 
7r/o at small momentum transfer and high incident en- 
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ergy. We denote these = Td( 7 tN —>► 7 rpN) and 
Tjj = T d (ttN —>► 7rfo N). The reactions are represented in 
FigsH](a) and (b). 

The up case has been studied at length. Its amplitude 
(c.f., Eq. ( 2 . 1 ) of Ref. jB[) is 

Tb = g P 7^K p (t 2 ) 2 1 is 13 e btl (Tnp ( 1 ) 

ra 2 - t 2 

where g p7T7T is the pnn coupling constant (g 2 /( 47 r) = 2.4), 
K p is the magnitude of the incident pion momentum in 
the p rest frame Q, b is the slope of the ttN elastic diffrac¬ 
tion peak, and a np is the tt p total cross-section. The 
invariants 813 , t\ and £2 are labeled in Fig. [T| (a). 

Similarly, the 7r/o production amplitude is 

~~ Q foTTTT o I - ^13^ 1( ^7Tp • ( 2 ) 

— t 2 

Choosing the average value of the /o tttt width of 
60 MeV, we obtain a numerical value ^/ 07r7r — 1.45 GeV. 
The other factors in Eq. 0 , relative to Fig. |T] (b) have 
the same meaning as in Eq. O- 

The p7T Deck background has been well studied, where, 
by background, we mean the amplitude before any final 
state interactions are included. We refer to Ref. [5| and 
extract what is useful in the present analysis. We work 
in the final pir (/o^r) center of mass frame; M is the 
invariant mass of this system. In the limit of forward 
production ( t\ 0) and large 8, Eq. dTJ produces the pir 
system predominately in an S-wave [4J, used in previous 
calculations (e.g., Ref. i). 

However, the J p = 1 + /o7r system is in an orbital P- 
wave. To address /q7t, we must extend the partial wave 
extraction calculations to finite values of t\ and 8. We 
present the complete calculation of these amplitudes else¬ 
where The important feature is that the higher par¬ 
tial wave amplitudes are of order t\/M 2 or M 2 /s with 
respect to the dominant S-wave. An immediate conse¬ 
quence is that /o7r P-wave production should have a no¬ 
ticeably smaller rate than the pir S-wave process, as is 
borne out in the complete calculation and exhibited by 
the COMPASS data, where the intensity of the /q7t peak 
at 1.42 GeV, is lower than that of the pir peak at 1.26 GeV 
by a factor of the order of a few 10 -3 . 

In the COMPASS experiment, the value of the square 
of the invariant total energy is s = 380 GeV 2 while 
the momentum-transfer t\ in the smallest bin is t\ G 
[—0.1, —0.13] GeV 2 . We are interested in values of M ~ 1 
to 2 GeV. Since \ti\/M 2 M 2 /8, the only relevant 

kinematic corrections come from the momentum trans¬ 
fer dependence. We choose to work at the fixed value 
t\ = —0.1 GeV 2 , and we checked that within the first 
t-bin (£1 G [-0.1,-0.13] GeV 2 ), our results do not vary 
appreciably. A convenient dimensionless expansion pa¬ 
rameter is 


(M 2 — m 2 ) ’ 


The J p = 1+ S-wave pit background amplitude is, to 
first order in @ 1 , 
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where E n and E p are the pion and p energies in the pir 
rest frame, and where y = p^/E^ is the pir phase space 
factor, Pn being the pion momentum in the pir rest frame. 

The J p = 1+ P-wave /o7r amplitude is, at the same 
order in @ 1 , 
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where E n , Ef 0 are the pion and /o energies, p n the pion 
momentum in the fair rest frame and, as above, y = 
Pn / E n . 

Equation (j5j) is a major clue to our investiga¬ 
tion. The right hand side contains the factor 
(3 M 2 + m 2 )/(M 2 — m 2 ) — Ef 0 /E n . This factor is neg¬ 
ative at low values of M (since m^/m^ > 3), but it 
vanishes near M ~ 1.38 GeV and becomes positive after¬ 
ward. Furthermore, if we give this term some small imag¬ 
inary part, its phase will switch suddenly from —180° to 
zero. This sudden and rapid phase variation is not a 
dynamical effect in the sense of a resonant phase, but 
it originates in the structure of the dynamical process 
by which the /o7r state is produced. Another interesting 
qualitative feature of Eq. (J5j) is that it grows in the region 
of interest (M ~ 1.2 to 1.4 GeV) and therefore tends to 
push a resonance peak upward in M. 

Keeping in mind the parameters introduced in Eqs. 0 
and 0 , our two J PC = 1 ++ amplitudes are 


f Tn ec k(p7t) \ _ 2iV2sN f 9p7T7T Epd^pTp^ \ 
v T D eck(f ott) ) (M 2 - m 2 ) V 9fo7T7T cr^p T fon ) ’ 
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where, T piv and X/ 07r can be read off from Eqs. (|2|) and 
([5|). The structure remains the same after we unitarize. 
The normalization factor N is irrelevant for present 
purposes and is taken equal to 1 here. 


Unitarization. For theoretical and technical details 
about multi-channel final state unitarization, we refer to 
the literature, in particular to Ref. 9| where the general 
analysis may be found and to Ref. [5| where a specific ap¬ 
plication is made. We recall that if S is the (two-channel) 
strong interaction S matrix, then the unitarized Deck 
amplitude 7x>, which we can write as a two dimensional 
vector as in Eq.(j6]), has a right hand unitarity cut along 
which it satisfies the relation 

T v + = ST V ~ 


©1 = 
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Tv + and Tv~ being the values of the unitarized Deck 
amplitude above and below the cut. 

Our basic assumption is that there is a single a\ res¬ 
onance whose (unique) second-sheet pole parameters we 
determine. Since we are dealing with a two-channel case, 
we parametrize the coupled pir and fair final state in¬ 
teractions (or rescattering) via this resonance. In order 
to do this, we introduce a K matrix, as in Eq.(3.14) of 
Ref. 


Analysis and results. Some salient points can be made 
short of a detailed fit to data. We first select appropriate 
values of the a\ parameters that provide a good global 
description. The COMPASS results fix these parame¬ 
ters more stringently than when we dealt only with the 
S-wave pir system (and other S-wave channels, such as 
K*K). Here, the acceptable mass and width of the ai, 
defined by the position of the second sheet pole , turn out 
to be quite restricted. Our analysis indicates that: 


( 3 ? 9192 \ 

si ~T S1 T I • ( 8 ) 

si—M 2 si-M 2 / 

The crucial tool to treat coupled channel final state in¬ 
teractions is a D-matrix, related directly to the S matrix. 
It is presented explicitly in Eq.(3.15) of 


D(M 2 ) 


1 ( 91 -# 2 (si - M 2 - a 2 C 2 ) \ 

Vo(M 2 )\g 2 gi(s 1 -M 2 -a 2 C 1 ) J 


(9) 


where a 2 = Q 2 + g 2 , C\ and C 2 are Chew-Mandelstam 
functions |Rj], and the energy denominator function 
V 0 (M 2 ) is 

Vo (M 2 ) = ( Sl - M 2 - g 2 C\(M 2 ) - g 2 2 C 2 {M 2 )) . (10) 


The function Vo (M 2 ) contains all the information (that 
we put in) on the coupled-channel pir-foTr strong inter¬ 
action. It is an analytic function which possesses the pir 
and /o 7 r branch cuts from [mp+m ^] 2 to infinity and from 
[m f o Am^] 2 to infinity. Its second sheet pole determines 
the nominal position and width of the a\ resonance. 

As in [5], the unitarized Deck amplitude with resonant 
rescattering corrections taken into account is 
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Here Tfi(M 2 ) is a two-dimensional vector and Tp(M 2 ) 
is the “background” Deck amplitude discussed above. 


Direct production contribution. In addition to its man¬ 
ifestation through final state interactions, the a\ may also 
be produced directly in a diffractive process up -A a\p. 
For direct production, we choose 

< i2 » 

where G represents the diffractive coupling to the i r to 
the ai, and f\,f 2 are the couplings of the a\ in the pir 
and foTr channels respectively. Our final amplitude is 

T(M 2 ) = Tp(M 2 ) + T d i r (s, M 2 ) . (13) 


In the one-channel case exemplified by p photo¬ 
production [ll|, the interference of these two terms can 
shift the apparent peak position of the p. 


M(a 1 ) - 1.40 ± 0.02 GeV, 

T(ai) ~ 0.30 ±0.05 GeV. 

These values of mass and width correspond to values of 
the parameters si ~ 2.002 GeV 2 and g\ ~ 0.732 GeV. 

The interesting parameter to vary is the ratio 7 = 
g 2 /gi in order to find the range of values that produce 
two peaks with appropriate characteristics: the /o 7 r peak 
occurs at higher mass than the pn peak, and the ratio of 
maximal intensities of these peaks, i.e. p 7 r//o 7 r, falls be¬ 
tween 1,000 and 500, as indicated by the available data. 
These requirements lead to negative values of 7 in the 
range [-0.1,-0.055]. In other words, the a\ couplings 
to pir and to /o 7 r have opposite signs. We choose as our 
central value 7 = g 2 /gi = —0.08. 



FIG. 2: The pn differential cross section as a function of the 
mass M, in three cases: background Deck, unitarized Deck as 
in Eq. CD- and the final result including direct production. 


The next step is to determine the amount of direct 
production necessary to fix the two peaks, one in p 7 r, 
the other in /q 7 t, at their desired positions, i.e. M = 
1.26 GeV for pir and M = 1.42 GeV for /q 7 t. Values such 
as Ga np fi = 120 and Ga np f 2 = 5.5 in Eq. (fl2|) ensure 
good positions for the two peaks, and this situation is 
stable when one varies the parameter 7 . The ratio of 
direct production to the background Deck amplitude is 
consistent with the value we obtained previously 0 in 
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our analysis of data at much lower energies, with smaller 
statistics. 



FIG. 3: The /o7r differential cross section as a function of the 
mass M, in three cases: background Deck, unitarized Deck as 
in Ea. (fill) , and the final result including direct production. 


It is interesting to see how the two peaks are built up. 
We plot in Fig. (|2]) the shapes of the pn intensity as a 
function of the energy M, for the various terms in the cal¬ 
culation. The pure Deck background does not produce 
a resonant shape. The unitarized amplitude shows effec¬ 
tively the ~ cos£ zero that appears in the one-channel 
case [5] (to which this problem is actually very close). 
Finally, direct production produces the observed peak, 
at the right position. A similar set of curves for the fair 
intensity is shown in Fig. (j3j). Notice that the form of 
the Deck background by itself appears to simulate a nar¬ 
row resonance peak at threshold (of course without any 
accompanying phase). 



FIG. 4: Three phase-differences between fon and p7r, in the 
1.2 to 1.6 GeV region. The middle curve corresponds to our 
“central” solution 7 = —0.08; the two others (marked) to 
other values in the range of interest. In all three cases, the 
other results of the calculation remain practically unchanged: 
peak positions and intensities. 


The separation of the positions of the two peaks is 


evident. The width of ai(1260) is about twice the width 
of ai(1420) in the calculation. The ai(1420) peak is also 
more symmetrical, with width about 0.14 GeV. The lower 
end of the fair intensity exhibits a (tiny) peak at around 
1.2 GeV owing to the zero emphasized in Eq. ([5]). 

We display in Fig. a set of phase differences be¬ 
tween the /o7r and pir amplitudes for three values of our 
parameter 7. We adopt 7 = —0.08 as our “central” value 
(subject to more refined analyses). We call attention to 
the rapid rise of over 100° in the phase difference in the 
mass region 1.2 to 1.3 GeV. A quantitative fit to the data 
would yield a precise value of 7 and of the branching ra¬ 
tio of the a\ into /o7r of the order of 10 -3 relative to the 
dominant decay mode a\ —>• pir. 

Summary. We find that the main features of the COM¬ 
PASS data, two mass peaks separated by ^ 160 MeV 
with significant relative phase motion, are fully compat¬ 
ible with a single a\ resonance. A detailed quantitative 
fit of the data with this formalism would lead to a new 
determination of the mass and width of the a\ and of its 
branching fraction into /q7t. 
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